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, , , Abstract 

X 

' The inclusive decay rate into pions of the charmed Ds meson is surprisingly 

larger than estimates expected from the W annihilation, adopting commonly 
used values of current-algebra up and down quark masses. We then go beyond 
this tree diagram and consider possible QCD effects that might cause such 
a large rate. There are two; the first one is related to the spectator decay 
cs ^ ss + ud followed by ss ^ dd , uu via two-gluon exchange box diagram. 
The second one is a gluon emission in weak annihilation for which the usual 
helicity suppression is vitiated: Dg W + g followed by — > ud, g — > dd, uu. 

These two contributions, however, turn out to be insufficient to explain 
data, implying that the puzzle could be understood if the up, down quarks 
have higher mass values. 

Furthermore, on the basis of experimental informations on the spectral func- 
tion P3tt{Q'^) deduced from the exclusive Ds — > Stt mode, the QCD sum rules 
also point to a higher mass for light quarks. 



The inclusive decay rate of the charmed Dg meson into non-strange ordinary 
hadrons ( mainly pions) is surprisingly large Its amplitude, being governed by the 
W annihilation mechanism at the electroweak tree diagram and directly related to the 
divergence of the axial current, is expected to be strongly suppressed by the familiar 
helicity argument ( similar to the suppression of vr — ez/ compared to vr ^ fiu), 
and/or by the partial conservation of the axial current (PCAC). Experiments do not 
confirm this expectation, however. 

We observe first that the Dg decay into pions cannot be described either by the 
dominant spectator mechanism ( both color-favoured and color-suppressed) or by the 
small penguin diagram, because the spectator constituent s of the Df is absent in the 
decay products. Therefore only the W annihilation mechanism can give rise - at the 
tree level - to the decays of into pions. The gluonic effects will be also considered. 
Experimentally, the inclusive branching ratio for Dg decays into non-strange hadrons 
( denoted in the following by Xud) can be estimated to be at least (1.65 ± 0.35)% 
from the Particle Data Group (PDG) For this lower bound, we only retain the 
three and five charged pion modes. All other modes with four, six, seven pions are 
disregarded because they come mainly from the quasi two-body modes Dg rj + Ti{p) 
and Dg ^ rj' + 7i{p) followed by the subsequent decays of rj and f]' into pions. Since r] 
and T]' have a large ss component, these modes must be attributed to the dominant 
spectator diagram and not to the W annihilation in which we are interested here. 
Let us remark that with the W annihilation, the decays of Dg into rj + tt and rj' + tc 
vanish by the conserved vector current (CVC) in the standard factorization approach. 
Another remarkable feature of the experimental data is the important fraction of Dg 
decays into non-resonant vr+vr"??"'" state. Its rate, which is already around 1/3 of the 
dominant spectator 07r+ one, is really intriguing. 

In the rest of the paper, we use two different methods - appropriated to inclusive 
and exclusive decays - to estimate the up, down quark masses. Both approaches 
converge to mass values higher than the ones estimated in literature. 



1 Inclusive Decay of Dg into pions: 

Compared to the pure leptonic rate Dg yU + z^, the inclusive decay of Dg into 
non-strange hadrons Xud is given, at the effective electroweak tree level, by: 

_ r(£'+ Xud) _ r(£'+ ud) _ ^ ^ J{xu,Xd) ml + ml 
""^^ - r(D+ ^ ^+z.) T{Dt^f.^u) (l-x,)2 ml ■ 

By considering the ratio Rd^ in Eq.(l), the large uncertainty in the decay constant 
fo^ can be avoided. The coefficient 3 comes from color. Here ai is the Bauer- 
Stech-Wirbel (BSW) phenomenological parameter taken from the QCD corrected 
effective Lagrangian first calculated by Gaillard-Lee , Altarelli-Maiani following 
the Wilson operator product expansion method 

G 

^eff = -J=V*gVud[ci{s-f^Lc){u-f^Ld)+C2{s-f^Ld){u-ff,Lc)] (2) 
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In the vacuum insertion approximation, known as the factorization method a la BSW 
p|, the relevant effective Lagrangian appropriated to our case can be written as 

G 

jO^eff = -^V*^Vudai{s^t^LC)H{u^/iLd)H (3) 

with ai = ci + -^€2 and the subscript H stands for color singlet hadronic currents 
i- 

From the general fit of charm non-leptonic decays, the value of ai = 1.26 is 
commonly used 

Finally J{xu, Xd) in Eq.(l) is the phase-space correction factor due to the quark 
masses ^ 

J{Xu,Xd) = [1 - ^^"^ "^"^ ] X{l,Xu,Xd) (4) 
•^d ~r Xy_ 

2 

"if, 



with A(l,x,2/) = ^{1 - X - yy - Axy, = where m„, m^, and M 

are respectively masses of up, down quarks, muon and charmed meson. 

We remark the similarity of the ratio -R^i^ with the hadrons/lepton ratio Rr in the 
tau's lepton decay 

Rr = T^T— r _ , . = 3\Vud\ F{xu,Xd) (5) 

1 (r+ —>■ Ut-c+u) 

with the phase space factor F{xu, Xd) given by 

F{x,y) = X{l,x,y)[l-7{x + y)-{x'^ + y'^)-Q{x-y)^ + {x + y){x-yy 

- Gxyix + y)] + I2[x'{l - 1/') | ^ I ^ ! I' ^] + ^ v)] (6) 

J(0,0) = F(0,0) = 1. 

Eqs.(l) and (5) are only exact at the tree level, i.e. when QCD effects - at the 
final up, down vertex- are not taken into account. When gluonic effects ( up to three 
loops) at the final quarks vertex are kept |^, then both Rr>^ and Rr in Eqs.(l) and (5) 
respectively have to be multiplied by a common correction factor[| G given by 

G = 1 + — Xd) + { — y^iiXu, Xd) + { — )^63{Xu, Xd) (7) 

vr vr vr 



* Let us emphasize that there is no confusion ( or double counting) possible between the QCD 
correction in the effective Lagrangian Eq.(2) symbolized by the ci^2 coefficients on the one hand, 
and the QCD corrected coefficient G at the u, d vertex on the other hand. 

The first one ci,2 are quantities that issue from the renormalization group equation that sums up 
large logarithmic enhancement due to gluons crossing the W line, i.e. gluons connecting the initial 
sc to the final ud. This is why r decay is not concerned with this effective Lagrangian ( gluons 
ignore r lepton). 

'I' The QCD corrected coefficient G in Eq.(7), on the other hand, operates only at the final state 
u, d vertex, therefore G is common to both Rr and Rd^ ■ 
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with 

(5i(0, 0) = 1 , 52(0, 0) = 5.2 , 53(0, 0) = 26.36 (8) 

In Eq.(7) non-perturbative QCD contributions |§] turn out to be tiny and are conse- 
quently neglected. 

For non-zero arguments Xu,Xd, the one-loop function 6i{x,y) has also been com- 
puted with a rather surprising feature Si{x,y) > 5i(0,0) recently confirmed [|T0|] . 
To our knowledge, 62{x,y) and 6s{x,y) for non-zero x,y are not yet computed. 

Now the first question that arises is which mass - current algebra mass or con- 
stituent mass - must be used for the up, down quarks? 

We argue for the first one ( current mass) due both to theoretical and phenomeno- 
logical reasons. Theoretically, in QCD, the divergence of the observable axial weak 
current is given in terms of current mass, as extensively discussed in Ref.[ll]. The 
constituent mass is only appropriate for the bound state problem not considered here. 
Let us also remark that the hadronization of u, d quarks into pions is conceptually 
different from the boundstate problem. 

At the phenomenological level, the use of constituent mass ~ 300 MeV would 
firstly spoil the excellent arguments in favour of QCD tests in r decays [0, |], |12 



indeed, the enhancement of about 20% through the QCD factor G in Eq.(7) would 
be diminished by a similar amount if constituent mass is used because of the phase 
space factor y) of Eq.(6). Secondly, for constituent mass, the ratio Rd^ of Eq.(l): 
would be in strong disagreement with data by two orders of magnitude. 

Once the principle of current algebra masses is retained, let us return to Eq.(l) 
and ask ourself which numerical value must be used for iriu and m^. 

This is because, contrary to in Eq.(5) which is insensitive to exact values for 
current masses, the ratio as given by Eq.(l) is very sensitive to them. 



Light quark masses have a reputation of being " not well measurable" [0, |T3l and 
even their indirect experimental determinations have so far not even been attempted 
111 , while there exists a huge number of theoretical estimates [0 ranging between 



4 MeV and 12 MeV for the average mass m of the m, d quarks. 



m = \{mu + md) (9) 



Since consensus holds [|T3l for the ratio r = "^'^ , ™'" = 0.28 ± 0.03, the problem resides 
only on m, more precisely on the running mass m(so) at the scale Sq of few GeV^. 

Returning to the left-hand side of Eq.(l), we first estimate the experimental value 
of Rd^ from the PDG [|l| to be 2.8 ±1.2 where the branching ration Br{Df Xud) = 
(1.65 ± 0.35)10-2 has been used together with Br{Df fi+u) = (0.59 ± 0.22) lO'^. 

Within the one standard deviation limit of data for Rds, Eqs.(l), (7) can only be 
satisfied with m around 38 MeV - a huge value much larger than those commonly 
expected [|T^. In principle, there is nothing wrong with such a large m value, although 
this would imply a considerably lower value for the — < tpip > condensate than the 
standard chiral perturbation theory could support. Therefore before taking seriously 
this crude 38 MeV value, we must ask ourself the next question, what could be the 
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other sources that may contribute to the substantially observed inclusive branching 
ratio Br{Ds X^d) ? 

There exists at least two possibilities. The first one is via the so-called Zweig 
forbidden rule as depicted in Fig.l: the dominant color-favoured decay mode c 
s + iud) when combined with the spectator s could induce, through a two-gluon box 
diagram, the chain cs — * (ss) + {ud) (qq) + {ud) where q stands for u and d quarks. 
It is the final state interaction at the quark level (appropriated to inclusive processes), 
in which the final state ss turns into qq. 

The second one, as depicted in Fig. 2, is a weak annihilation accompanied by a 
gluon emitted from the initial s and c quarks bound inside Dg : mechanism, first 



proposed in Ref. fl^ to vitiate the helicity suppression and recently being reexamined 
in details |]T5[, could a priori yield a large Dg — > X^d rate ( since the width is no 
longer suppressed by ^). 

Our task now is to compute these two contributions of Fig.l and Fig. 2, that we 
called respectively final state strong transition (FT) and gluonic weak annihilation 
(GA). 

1) For the first one (FT), the decay rate can be written in the form 

Tft{cs {qq) + {ud)) = Tsp{c -> s + {ud))P{ss gg ^ qq) (10) 

where Tsp{c s + {ud)) is the familiar spectator inclusive rate and P{ss gg ^ qq) 
is the transition probability for the ss pair in the final state transforming into a qq 
pair through a box diagram. 

The first term in Eq.(|lOD rsp{c ^ sud) is given by: 

n'^m^ In SI 

Tsp{c ^ s + {ud)) = 3a?— -^|K.Kdr/(x„ x., x,)[l - -^(tt^ - -)] (11) 

In Eq.(|ll]), the phase space factor I{x,y,z) is known with /(x, 0, 0) = 1 — 8a; -|- 
8x^ — — 12x'^lnx. The QCD correction factor C = — ^(tt^ — ^) corresponds only 
to massless s, u, d quarks and can be decomposed into two parts, the "upper vertex" 
associated to the cs part is — ^(vr^ — ^) and the "lower vertex" associated to ud part 
is For massive s, u, d quarks, both "upper" and "lower" parts have already 

been computed |^ , the explicit analytic expression for the "upper" part is also given 
rD|,|T§. 



We now compute the box diagrams ( there are two, with crossed gluons). The 
dimensionless box diagram amplitude corresponding to s{pi) + s{p2) — > q{p3) + q{p4) 
can be conveniently written as: 

^ 3 Q 

with 

Q' = (Pi+P2)', t={pr-pzf (13) 
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The dimensionless quantity B{Q^, t, m^) coming from loop integration has the follow- 
ing representation ( we neglect ) 

B{Q ,t,m,) = l dyj^ dzz{l - z) D^^.^y, ,,Q2^t,m^^) 

(14) 

where 

D(x, y, z, g^ t, m^) = -Q'x(l - x)(l - zf - ty(l - y)z'' + mlyz'' (15) 
Its explicit expression is given by: 

B{Q\t,ml) = ln\-'^) + —^C{Q\t,ml) 



g2 + 2i-4m2^ + 2^-2™' 



with 



where is the Spence or dilogarithmic function. 

Once the box amplitude A for ss — > qq is known, the transition probability Y = 
P{ss gg ^ qq) in Eq.(lO) is then obtained by: 

Since / ^ / ^5'^{Q —ps — P4) = |, the quantity Y can be conveniently rewritten as 

Y=(^ylA(Q\ml) (18) 

The range of variation for Q^: 

4m^<Q'<— (m, + me)' (19) 
rric 



is deduced from < P < {rric — m^)^ where / is the invariant mass of the ud pair 
issuing from the spectator mechanism c ^ s + (ud) . 

The angular integration over the variable t in Eq.(^) is done numerically, and 
AiQ^,m1) turns out to be inside 0.36 < A{Q^,m'l) < 0.95 for inside Eq.(|l|), 
using nis = 0.15 GeV, rric = 1.45 GeV. 

With these values put into Eqs. ([TO|) , ( pT]) and (0), the FT contribution yields a 
branching ratio Br{Ds pions )ft ^ 0.68% which is still far from the experimental 
data, which at least equals to (1.65 ± 0.35)%. 

2) For the gluonic weak annihilation (GA) of Fig. 2, the amplitude can be directly 
taken from Eq. (4) of Ref . [|I4[ with the following substitution: 

q^h + h , gse){q) ^ glu{k{)Y^v{h)j^^-^^, (20) 

where ki and k2 are the q and q momentum, and j denotes the color index. 
We have 

^ = ^ ^ j^^y u{qi)Y{^ - l5)v{q2)u{ki)-f -v{k2), (21) 

with ^ 

R^uiP, q) = J—^{FA{q^,Pu - {q.p)g,,u) + iFye^.^apP^'q^}, (22) 

where Fv,a are the two dimensionless form factors somehow reflected the wave func- 
tion at the origin of the cs bound state, which in turn is proportional to the decay 
constant f^^ via the Van Royen- Weisskopf formula. The tensor R^^ is reminiscent 
of the standard fermionic triangular loop describing 'yn'^ or it" 77, taken 

as an example. 

Here p = gi + g2 + ^1 + ^2 with = M^. Neglecting m^^ compared to M^, 
the four-body phase space integration of ( spin and color summing up) can be 
simplified, and the decay rate is computed to be: 

ToA{Ds^ud + qq) = a f ' , L 3 / / dt^t,{M' - t^f 

IZTT'^M Jo 



2) 



with ti = (gi + ga)^ , ^2 = (^1 + ^2)^ 



The double integration over ^1,^2 can be done analytically and the result is 



2G^\V:M^al{F^ + Fl)M\n^ , 10,^^ 71. 



TcAiD, ^ud + qq) = af ' ^ ^^^3^ {- + -ln2 - -). (24) 



Following Ref. ||14||, we will take 
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Numerically, it turns out that the Tqa rate is extremely small, the corresponding 
branching ratio Br{Ds pions) via gluonic weak annihilation mechanism is at most 
(0.1)%. In Eqs.(12) and (24) we take a, = 0.3. 

Putting altogether both FT and GA contributions, the inclusive branching ratio 
Br{Ds pions) due to these gluonic processes cannot exceed 0.8% and is still far 
from the observed inclusive branching ratio larger than (1.65 ± 0.35)% as estimated 
from PDG §. 

The difference between these two numbers must be attributed to the pure W 
annihilation tree diagram I, from which at the one standard deviation lower bound 
of data, we get m ~ 22 MeV using again Eqs.(l) and (7). 

This value, although smaller than the crude 38 MeV obtained above ( when the 
gluonic backgrounds are neglected) is still at least twice as large as the common 
estimates. 



2 Exclusive three pion mode: Mass determination 
from QCD sum rules 

We now consider the exclusive decay mode —>■ tt+tt+tt^ which allows us to extract 
the spectral function p{Q^) that will be in turns exploited in the QCD sum rules to 
obtain m. 

The starting point is the QCD sum rule for the two-point correlator of the 



divergence of the axial- current, put in the form [11 



POO 

m\so) = H-\w, So) / dQMQ', So)p(Q') (26) 
Jo 



We also consider its finite energy version closely followed Ref. 

^V,^o) - (07^)^S^[l + ^2(-o) + 2c4<O4>/.^ r 

{flml+ rfQW)} (27) 
Jo 

In Eq.([26|), w(Q^, sq) denotes weight function and H{w,so) is defined by the large 
behavior of the two-point correlator, it has a perturbative QCD part and a non- 
perturbative part parametrized in terms of vacuum condensates. Also the two-loop 
expression for R2{so) as well as the dimension-4 condensate C4 < O4 > in Eq.(0) can 
be found in Refs.[jl|, H. 

These types of sum rules have been extensively used in the literature to estimate 
quark mass and as explained in Ref. , the main problem in these estimates is not the 
sum-rule technique itself but rather the complete absence of experimental information 
on the magnitude of p{Q^) beyond the one-pion contribution. The later f^m^ is 

^ The W annihilation and the two gluonic contributions have different decay products, therefore 
they contribute incoherently to the Dg pions rate, without interference. 
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singled out in Eq.(^) with /^r ^ 132 MeV. Keeping only this pion contribution, one 
finds ||Tl|, |18[ a lower bound m(l GeV) > 4(5) MeV due to the positivity of piQ"^). 



Fortunately, we show that the decay mode Svr provides precious information 

on p{Q^) at fixed = M^. Indeed, the decay amplitude tt+tt+tt" is given by 



G_ 

where the most general matrix element < 37r|y4^|0 > can be expressed in terms of 



three form factors [19|, the same occurred in r — > u + Sn. However, for Dg —> Stt, only 



one dimensionless form factor pO) associated to the divergence of the axial current is 



involved and will be denoted by F{si, S2) 0: 

< 7r+(pi)7r+(p2)7r-(p3)|^l0 >= MF{s,,S2) (29) 
with Si = {Q - , i = 1, 2. 



In terms of -F(si, S2), the Stt rate is given by Ref.pO[] 



T{D, 3n) = a f'^^^f M'fl K, (30) 



where the dimensionless constant K is obtained by integrating the squared form factor 
\F{si, S2)P over the whole Dalitz domain: 

^^m^I I dsids2\Fis„S2)\' (31) 
Ami <si< (M - m^f 

Smin{Si) < S2 < Smaxi^i) (32) 



^min,max 



.-^i) ± 



2 V 2 

The numerical value of K that we can extract from Eq.(^) has two sources of errors, 
the first one is related to the experimental errors of the non-resonant branching ratio 
Br{D^ 7r+7r+7r^)ArR = (1.01 ± 0.35)10"^, the second one is the uncertainty of the 
decay constant /d,. Fixing fo, = 280 Mel/^, we get ^ K = 0.486 ± 0.168 where 
errors come from the ones of the experimental Dg — > tt+tt+tt^ branching ratio. Let 
§ fD, = 280 ± 70 MeV, then we have: 

0.27 ± 0.06 < K < 0.74 ±0.17 (33) 

Now the crucial point is that the constant K can be directly related to the spectral 
function ps^iQ'^) entering in the QCD sum rule Eqs.(26),(27). At = M^, we get: 

PsAQ' = M') = (34) 



§ Our dimensionless form factor F{si, S2) is related to the F4 form factor of Ref. |i9|| by F{si, S2) = 

MF4(S1,S2,Q' =M2) 
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Eq.(34) is obtained if we remind that the spectral function is defined by 

P(Q') = 7^E(2^)''^'W < ^I^M^IO > p (35) 

and in particular for p^.-niQ'^)'- 
P^^ = ll I I ###7A^(2vr)^5^(Q-pi-p.-P3)| < SvrlAlO > (36) 



27: J J J 2Ei 2E2 2E3 (27r)9 

In Eqs.(33) and (34), the constant K and consequently the p^^^^M'^) are only related 
to the charged — > tt+tt+tt" mode. The neutral mode D'^ —>■ 7r°7r°7r+ by isospin 
consideration ||2^ is presumably 1/4 of the charged one. Hence the total charged and 
neutral must be 5/4 of Eq. (^4]) . 

As shown in Eqs.(30) and (34), from Ds 37r data, we get information for the 
spectral function p{Q^) at only one fixed value of = M^, while the QCD sum rules 
Eqs.(26) or (27) need the whole range of in p{Q'^). Fortunately, we also know 
that at large Q^, perturbative QCD gives |]TB| explicit analytic expression for piQ"^) 
in terms of one unknown parameter (m)^: 

' 27r2(i/^g2/A2)iL 3 vr ^ ^ ^ 

Since the normahzation at Q"^ = for psj^iQ^ = M'^) is fixed by ^^M"^^ let us fol- 
low the usual procedure by adopting the simplest duality ansatz parametrization 
for p(Q2): 

3 7r 

which reduces to piJS/P') = as it should be. 

Putting now Eq.(38) into the sum rule Eq.(26), and taking for K the one standard 
deviation most conservative lower limit in Eq.(33), we then obtain for m(l GeV) a 
quite large value 20 MeV, consistent with our previous result from the inclusive decay 
Ds — >■ Xud- 

A different parametrization for p{Q'^) can be obtained by making use of the shape 
of the low energy end of the 37r-threshold p^if) calculated in a recent paper |^ using 



lowest order chiral perturbation theory. Following Ref . ||22|| , we will make use of the 
finite energy sum rule Eq.(27) with p(Q^) given by 

PiQ') = A^^P...(Q^) (39) 

where PhadiQ^) is taken from Eq.(32) of Ref.|2^ encoding the presence of two n' 
resonances, the 1300 MeV and 1770 MeV. 

Here A is an overall normalization which, as emphasized by the authors of Ref.[|ll|], 
can only be determined by experimental informations and not by theoretical calcula- 
tions. Matching Eq.(34) to Eq.(39) at = M^, one can fix A. It turns out that our 
overall normalization A provided by Ds —>■ Sir decay rate is much larger than the one 
given in Ref. p2[ , and that is the crucial reason for getting high value of rh. 
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3 Summary and Conclusion 



The starting point is our observation that the decay rate of the charmed Dg meson into 
pions is surprisingly larger than common estimates. The amplitude being governed 
by the W annihilation is usually expected to be negligible by the partial conservation 
of the axial current ( PCAC) at relatively high energy. 

We then go beyond the tree W annihilation diagram and consider two possible 
contributions. The first one, called finite state strong transition ( FT) coming from 
the dominant spectator decay cs — ss + ud followed by the Zweig violating rule 
ss qq {q = u, d) via the two-gluon box diagram. The second one, called gluonic 
weak annihilation ( GA), is a genuine weak annihilation accompanied by an emission 
of one gluon from the quarks c, s bound inside Dg, the mechanism that invalidates 
the helicity suppression. 

We find out that these gluonic effects together cannot explain the large inclusive 
Dg decay into pions, and consequently we suggest that W annihilation (tree diagram) 
also contributes substantially to the rates. This would imply that m, the averaged 
M, d mass, could be around 22 MeV ^ few times as large as the usual estimates ||l3i . 

One might argue that our quark parton type of analysis for inclusive decays of 
charm could only be accurate within a factor of two. However, by analyzing the 
exclusive Ds —>■ Sir mode, we obtain the normalized spectral function p3n{Q'^), which 
is in turns exploited in the QCD sum rules. Again, by this completely different 
method, a higher value than common estimates is obtained for rh. 

Our lower bound of 22 MeV for rh depends, on the one hand, on the experimental 
data of both inclusive and exclusive Dg decays into pions, and on the other hand, 
on the theoretical factorization method a la BSW. To obtain such a mass value, we 
carefully take into account experimental errors as well as the uncertainty of the decay 
constant fn^- 

Independent of our proposition for a high m as a solution to the substantial 
rates of Dg into pions, the understanding of the origin of such data constitutes a very 
interesting problem in its own right. The confirmation of experimental data is equally 
important, on the other side. 

Finally we would like to emphasize the similarity as well as the complementar- 
ity between r lepton and charmed Dg meson decaying into ( three and more) pions: 
While r probes the dominant spin 1 and the small spin axial spectral functions in 
the whole range of the momentum transfer < < ( to separate them is a 
challenging experimental problem Jill), probes directly and easily the spin 

spectral function at only one value = This latter is the key for the determi- 

nation of the light quark mass, a fundamental parameter " not well measured" in the 
standard model. Its importance in chiral symmetry breaking from both perturbative 
and nonperturbative aspects is well-known. 
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